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Single crystals of MgCNi3, with areas sized up to 1 mm 2 , were grown by the self flux method 
using a cubic anvil high pressure technique. In low applied fields, the dc magnetization exhibited 
a very narrow transition into the superconducting state, demonstrating good quality of the grown 
crystals. The first critical field H&, determined from a zero temperature extrapolation, is around 18 
mT. Using the tunnel - diode resonator technique, the London penetration depth was measured with 
no applied dc field and the Campbell penetration depth was measured with the external dc fields 
up to 9 T for two different sample orientations with respect to the direction of applied magnetic 
field. The absolute value of the London penetration depth, A(0) = 245 ± 10 nm was determined 
from the thermodynamic Rutgers formula. The superfluid density, p s = (A(0)/A(T)) 2 was found 
to follow the clean isotropic s-wave behavior predicted by the weak - coupling BCS theory in the 
whole temperature range. The low - temperature behavior of the London penetration depth fits the 
BCS analytic form as well and produces close to the weak - coupling value of A(0)/fcflT c = 1.71. 
The temperature dependence of the upper critical field, H C 2, was found to be isotropic with a slope 
at T c of -2.63 T/K and H C 2(0) ~ 12.3 T at zero temperature. The Campbell penetration depth 
probes the vortex lattice response in the mixed state and is sensitive to the details of the pinning 
potential. For MgCNi3, an irreversible feature has been observed in the TDR response when the 
sample is field-cooled and warmed versus zero-field-cooled and warmed. This feature possesses a 
non-monotonic field dependence and has commonly been referred to as the peak effect and is most 
likely related to a field - dependent non - parabolic pinning potential. 

PACS numbers: 74.70.Dd, 74.25.N-, 74.20.Rp, 74.25.Wx, 74.25.Op 



The announcement of superconductivity in the in- 
termetallic compound MgCNi3 has generated a great 
amount of excitement since its discovery in 200 1.^ This 
material has gained so much interest because it is a su- 
perconductor with a transition temperature near 7 K 
and it shares the same perovskite structure as that of 
the high-T c cuprates but with the O atoms replaced by 
Ni. After the realization of these facts, many began to 
consider the possibility that this material could bridge 
the gap between conventional superconductivity in inter- 
metallic compounds and unconventional superconductiv- 
ity in high-T c superconducting oxide perovskites. This 
material has also been described as being a three di- 
mensional analog to the two dimensional family of boro- 
carbide superconductors. A possible scenario, brought 
about by predictions made for high-T c superconductors, 
is that the superconducting state arises due to interac- 
tions involving ferromagnetic spin fluctuations from the 
large Ni concentration. Band structure calculations indi- 
cate that the abundance of Ni in this compound places 
it near a ferromagnetic instability^^ and the existence 
of such a peak is confirmed by both photoemission and 
x-ray absorption experiments.^ 

Although some have made claims of observations of un- 
conventional superconducting properties, there are many 
measurements which indicate usual BCS behavior. C 13 



NMR investigations have found that the nuclear spin- 
lattice relaxation rate l/ 13 7i exhibits the typical expo- 
nential behavior expected for s-wave superconductivity.^ 
One point contact spectroscopy study concludes conven- 
tional weak - coupling BCS s-wave superconductivity,^ 
while another tunneling junction measurements have in- 
ferred the magnitude of the superconducting gap from an 
observed zero bias conductance peak and have found it 
to be larger than that of the weak coupling BCS value, 
which has led them to make the conclusion that the elec- 
tron coupling in this material is strong!^ A carbon iso- 
tope effect has been observed in this compound, indicat- 
ing that the carbon-based phonons do play an impor- 
tant role in the superconductivity and is thus in sup- 
port of the conventional BCS phonon mediated model 
of superconductivityP Muon spin rotation studies have 
found evidence for BCS behavior in the superconduct- 
ing gapP Measurements of the specific heat are indica- 
tive of a fully gapped superconducting state but they 
do not seem to be in agreement on the s trengh of the 
coupling or the effects of spin fluctuations! 10 * 11 ! Electri- 
cal transport measurements have found that the normal 
state resistivity follows a conventional electron-phonon 
scattering model and that H c2 near T c is linear. They 
have used these findings to conclude that MgCNi3 is 
a conventional BCS superconductor!^ Previous tunnel 
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FIG. 1. (Color online) Optical microscope images of MgCNi3 single crystals. An as grown melted lump with a mixture of 
single-crystalline MgCNi3 and some fluxes are shown in the left upper corner. After crushing the lump, a large number of 
crystals with the sizes up to 1 mm 2 were found. Upper right frame shows the hkO reciprocal space section determined by XRD 
of single crystal MgC0.92Ni2.88- 



diode resonator experiments on this material have been 
performed on polycrystalline samples and powders and 
it was found that the low temperature behavior was 
quadratic ,121 which would point to the existence of nodes 
in the superconducting gap function. However, it may 
have been the case that inter-grain interactions in those 
samples studied may have influenced the data. Most of 
the scattered physical properties and theoretical calcula- 
tions for MgCNi3 are reviewed in Ref.P^ L ater, with the 
appearance of the first single crystal datapEStHI jt be- 
came clear that there were some contradictions regarding 
the physical property measurements obtained on poly- 
crystalline and single crystalline samples. In addition, 
the recently observed peak effect and dynamics of vortex 
matter in MgCNijf^ requires further detailed investiga- 
tions. In this sense, MgCNi3 single crystals, which obey 
a simple perovskite cubic crystal structure, provide an 
interesting possibility for further magnetic studies. 

In this paper, high-pressure crystal growth and pre- 
cision measurements of the magnetic penetration depth 
on bar shaped MgCNi3 crystals for two different sample 
orientations in fields from to 9 T are reported. The 
superfluid density is constructed from the zero field pen- 
etration depth and this data has been shown to agree well 
with the isotropic BCS s-wave superfluid density model. 
The if c2 (T) curve is also constructed for fields applied 
in two different directions and it is found to be linear 
near T c and also isotropic, indicating that the change in 



T c due to demagnetization effects from the sample shape 
are negligible. The penetration depth in field, consisting 
of London and Campbell components, shows a very in- 
teresting hysteresis when the sample is zero-field-cooled, 
field-warmed and then field-cooled, most likely due to a 
vortex lattice response referred to as the peak effect. 

Due to the high volatility of Mg and the relatively 
poor reactivity of C, it is extremely difficult to synthesize 
single phase samples of MgCNi3 , even in polycrystalline 
form. The synthesis of single crystals is not possible in 
an open system; however, it can be done under high pres- 
sure, as was first demonstrated by Lee et a/.P^ As both 
methods show, the superconductivity in this material is 
very sensitive to the details of heat treatment and final 
stoichiometry. Amos et a/P^ reported that different C 
contents in MgC^Nia polycrystals caused different cubic 
cell parameters: a increased from 3.795 to 3.812 A as 
x varied from 0.887 to 0.978. In addition, the supercon- 
ducting transition temperature T c sensitively depends on 
the real C content and decreases with increasing content. 
In contrast to polycrystalline MgCNi 3 , the Ni site was 
partly deficient in single crystals synthesized under high 
pressure conditions.^ 

Here we report our successful growth process for 
MgCNi3 single crystals together with their structural 
and superconducting properties. The single crystals of 
MgCNi 3 were grown at ETH Zurich using cubic anvil 
high-pressure and high-temperature techniques. The 
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FIG. 2. (Color online) Magnetic susceptibility x as a function 
of temperature T for various magnetic fields H. 



mixture of Mg, C, and Ni powders in a molar ratio 1:1:3 
were placed inside of a BN crucible with the inner diam- 
eter of 6.8 mm, and the length of 8.5 mm. The heating 
element is a graphite tube. Six anvils generate pressure 
on the whole assembly. In a typical run, a pressure of 3 
GPa is applied at room temperature. While keeping pres- 
sure constant, the temperature is ramped up within 2 h 
to the maximum value of 1600-1700 °C, and is kept sta- 
ble for 1 h and then slowly cooled to room temperature. 
The high pressure was maintained constant throughout 
the growth and was removed only after the end of the 
crystal growth process. The final product was a melted 
lump with a mixture of single-crystalline MgCNi 3 and 
some fluxes (see left upper corner image in Fig.[TJ. Af- 
ter crushing the lump, the single crystals with various 
shapes and of sizes up to 1 mm 2 were mechanically ex- 
tracted (Fig. [I). 

The quality of the crystals was checked by using a sin- 
gle crystal x-ray diffractometer equipped with a CCD 
area detector (Xcalibur PX, Oxford Diffraction), which 
allowed us to examine the whole reciprocal space (Ewald 
sphere) for the presence of other phases or crystallites 
with different orientations. As it is clearly seen in the up- 
per right frame of Fig. [l] no additional phases, impurities, 
or intergrowing crystals were detected by examination of 
the reconstructed reciprocal space. The crystal structure 
was determined by a direct method and refined on F 12 , 
employing the programs SHELXS-97 and SHELXL-97P 
All atomic positions were found by a direct method. Af- 
ter several refinement cycles the correct crystallographic 
composition was determined and the final R factor was 
1.8% indicating the high quality of the structural model. 
The occupation parameters for the Mg, C, and Ni were 
found to be 1: 0.92, and 2.88, respectively. Thus, ac- 
cording to the structural analysis, the more appropriate 
chemical formula for our crystals is MgC0.92Ni2.88- Single 
crystal analysis confirmed the cubic structure with lattice 
parameter a = 3.7913(1) A. This value of lattice constant 
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FIG. 3. (Color online) Analysis in order to extract the lower 
critical field H c \ from m(H) measurements of the studied 
MgCNi3 single crystal. Upper panel: As measured mag- 
netic moment m versus the applied magnetic field H. Lower 
panel: Magnetic induction BV versus the applied magnetic 
field along the afe-plane. The inset presents the determined 
H c i as a function of temperature. We estimate the zero tem- 
perature value at jj,oH c i(0) ~ 18 mT. 



a is slightly smaller than that of observed in MgCNi2.s 
(a = 3.812 A, T c = 6.7 K) single crystals grown at the 
pressure of 4.25 GPa and temperature 1200 °CP^Lee et 
al. also note that crystals grown under a pressure below 
3.5 GPa had C deficiencies.^ The present data confirm 
this observation. However, in our high pressure growth 
conditions, besides the C deficiency, the Ni site was also 
partially deficient and thus the resulting T c is reduced 
more. For various growth batches, T c varies between 6.4 
and 6.8 K. For measurements of superconducting proper- 
ties, clean, fiat single crystals with sizes of a few hundred 
micrometers were selected. 

The susceptibility x of a flat, plate-like single crystal 
with approximate dimensions 0.8 x 0.4 x 0.1 mm 3 was 
measured using a Quantum Design MPMS magnetome- 
ter as a function of temperature for various magnetic 
fields applied along the planar sample. Field-cooled (fc) 
and zero-field-cooled (zfc) temperature dependent mea- 
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surements in low fields are shown in Fig. [2] The transi- 
tion to the superconducting state in low fields is very nar- 
row, demonstrating the good quality of the single crys- 
tal. Supplemental magnetization curves were recorded 
for this single crystal in the temperature range between 
2 and 5 K. From such measurements, the lower criti- 
cal field H c i was determined using a procedure discussed 
elsewhere For this the magnetic induction B was deter- 
mined from the measured magnetic moment and plotted 
as a function of magnetic field (see Fig. [3]). Due to the 
uncertainty of the sample volume V, only the product 
BV was calculated and plotted according to 



B = Hq(M + H) = fi (m/V + H) 



(1) 



Since B — in the Meissner state it is possible to cal- 
culate from the data of m(H), the field above which this 
equality is invalid. The sudden increase from zero oc- 
curs due to the penetration of vortices. The resulting 
tM)H c i(0) ~ 18 mT is consistent with a magnetic penetra- 
tion depth of ~ 200 nm, invoking a n ~ 100. All MgCNi 3 
single crystals investigated in this work did not show any 
traces of ferromagnetism, in contrast to a recent report 
of ferromagnetic domains coexisting with superconduc- 
tivity in carbon deficient MgCNiaP^l 

The single crystalline MgCNi 3 sample was studied us- 
ing a tunnel diode resonator (TDR) circuit technique. A 
detailed description of the application of this technique to 
study London and Campbell penetration depths in super- 
conductors can be found in Ref . and The principle 
elements of the setup consist of an LC self-oscillating cir- 
cuit supported by a tunnel diode. The tunnel diode has a 
heavily doped and extremely thin (10 nm) p-n junction, 
which gives it useful properties not common to ordinary 
diodes. The IV curve contains a region of negative dif- 
ferential resistance and when the diode is biased to this 
region, it acts as an AC power source for the tank cir- 
cuit. The tank circuit oscillates with a natural resonance 
frequency of /o = 1/2tt^/LC, which is very near 14 MHz. 
The sample to be studied is mounted on a sapphire rod 
and inserted into the inductor coil of the oscillator. The 
sample changes the resonance frequency of the circuit 
through its interaction with the ac magnetic field of the 
coil, which is on the order of 0.1 /iT. This small value of 
the excitation field of the coil ensures that its effect on 
the state of the sample is negligible and hence this tech- 
nique is non-perturbative. For a superconductor below 
its critical temperature, the ac magnetic field of the coil 
has a characteristic decay length, commonly referred to 
as the London penetration depth A, which is a function 
of temperature. The measured change in frequency, A/, 
is proportional to the dynamic magnetic susceptibility of 
the sample. This susceptibility may be written in terms 
of this penetration depth and a characteristic radius of 
the sample R, which is calculated using a procedure given 
in reference!^ So we have, 



Af(T) = -G X (T) = G 







1 tanh 




R 1 





(2) 



where the geometry dependent calibration factor is ex- 
pressed as G ~ foV s /2V c (1 — N), V s is the sample vol- 
ume, V c is the effective coil volume and N is the demag- 
netization factor of the sample. The factor G can be 
measured directly by extracting the sample from the coil 
at the lowest temperature of the experiment. Since the 
effective radius of the sample, R, is much greater than the 
penetration depth, A, this expression can be rewritten so 
that changes in the resonant frequency are proportional 
to changes in the penetration depth 



A/ (T) cx AA (T) 



(3) 



The most valuable feature that this technique has to 
offer is not the ability to measure the actual value of the 
penetration depth, but rather its variation with tempera- 
ture to great precision, AA = A(T) — A(T min ), with T min 
being the minimum temperature that can be reached dur- 
ing the experiment. The noise level of the system used 
for this experiment is ~ 0.1 Hz/hour, which combined 
with the natural resonance frequency of the system of 
14 MHz corresponds to a resolution on the order of parts 
per billion. This level of precision allows for the measure- 
ment of AA to a single Angstrom. The circuit assembly 
is mounted inside of a 3 He refrigerator that is lowered 
into the bore of a superconducting solenoid allowing for 
the application of dc fields up to 9 T in addition to the 
extremely small ac field supplied by the TDR. 

The superfluid density, p s , is an important quantity 
that can be related to the gap structure of a s uper- 
conducting through the London penetration depthP^^ll 
If the zero-temperature value of the penetration depth, 
A(0), is known, then the superfluid density can be con- 
structed from AA as. 



Ps{T) = [mt)) = 



AA(T) 

T(oT 



(4) 



where AA(T) is the measured variation of the London 
penetration depth, measured by using a TDR and by 
applying the calibration procedure described previously. 
Without a direct measurement, the most reliable proce- 
dure to evaluate A(0) is to use thermodynamic Rutgers 
formula that can be written asj^ 



dp 



T^rTc 



$ 



AC 



dH c2 
dT 



-A 2 (0) 



(5) 



T-s-T c 



Taking the measured slope, \dH C 2/dT\T^T c = 
2.63 T/K and the jump of electronic specific heat at 
T c , AC = 129 mJ/(mol(Ni) K), 11 and using the itera- 
tive procedure descrived in Ref. [2BJ we obtain A(0) = 
245 ± 10 nm, which compares reasonably well with 
A(0) = 232 nm determined from muon spin rotation 
measurementsP The symbols in Fig. 4 show the data 
(with the symbol size representing the ±10 nm error) 
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FIG. 4. (Color online) The superfluid density, p 3 {T), con- 
structed from the London penetration depth measured by a 
tunnel diode resonator (red circles). Symbol size represents a 
±10nm error bar. Expectations for the single - gap BCS su- 
perconductors are shown for clean - limit s— wave (solid blue), 
clean - limit d— wave (dashed green), dirty - limit s— wave (dot 
blue) and dirty - limit d— wave (dashed - dot green). Inset 
shows low - temperature variation of the London penetration 
depth and a weak - coupling BCS isotropic s-wave fit with a 
fixed A(0) = 245 nm and the gap, A(0), as a free parameter. 



and the lines show curves expected for single gap BCS 
superconductors in the clean limit s— wave (solid blue), 
clean limit d— wave (dashed green), dirty limit s— wave 
(dot blue) and dirty limit d— wave (dashed - dot green). 
Clearly, the clean limit weak coupling s-wave curve de- 
scribes the experimental data almost perfectly in the full 
temperature range. The inset in Fig. 4 zooms into the 
low-temperature region showing exponential saturation 
of the superfluid density approaching T = 0. More- 
over, if we use the measured A(T) = A(0) + AA(T) 
and fit it to the low - temperature expansion, AA(T) = 
A(0)v/7rA(0)/2A: B Texp(-A(0)/fc B r), where A is the 
maximum gap value at T = 0, in the "low - temperature" 
range of T < T c /3 we obtain an almost weak - coupling 
value for the gap to T c ratio, A(0)/£; s T c = 1.71. The 
BCS weak - coupling value for isotropic s-wave supercon- 
ductor is 1.76. Altogether our results convincingly es- 
tablish MgCNi3 to be a weak - coupling isotropic s-wave 
superconductor. 

Next we discuss the measurements of Campbell pene- 
tration depth in finite applied DC magnetic field. Tem- 
perature sweeps done in applied fields up to 9 T, where se- 
lected curves can be seen in Fig. 5, allow for the construc- 
tion of the H — T phase diagram, which is shown for two 
different sample orientations with respect to the applied 
magnetic field. The single crystal sample of MgCNi3 that 



was studied was a rectangular bar having approximate 
dimensions of 0.40 x 0.43 x 0.73 mm 3 . The two sample 
orientations about which both the ac and dc magnetic 
fields were applied are parallel to the long sample axis and 
perpendicular to the long axis with the fields being also 
along one of the principle axes of the sample. Fig. 5 shows 
that the H c2 curves for the single crystalline MgCNi 3 are 
isotropic. By analyzing these results within the Helfand 
and Werthamer theory we obtain a slope at T c of -2.63 
T/K and H c2 (0) « 12.3 T. The values obtained using this 
analysis are in excellent agreement with those obtained 
by another group performing resistivity measurements on 
single crystalline MgCNi3 in applied fields.^ 

Selected runs of TDR frequency shifts vs. tempera- 
ture performed in various applied fields and converted 
into susceptibility are shown in Fig. 6. For each run, 
the sample was cooled in a low field, with the first run 
cooled in zero field, and then the target field was applied 
after the sample had been cooled to the base tempera- 
ture. The resulting curve is independent of whether or 
not the sample is cooled in zero field or the previous 
field run value. The sample was then field- warmed and 
field-cooled twice. Notice from Fig. 6 that the initial 
zero-field cooled and field-warmed portions of the curve 
are irreversible, denoted by Xirr- This irreversibility is 
believed to be related to a response of the vortex lattice 
to the applied magnetic field and is a signature of the 
non - parabolic pinning potential.^ 

It should be noted that when considering the mag- 
netic penetration depth of a superconductor in applied 
dc fields, there are two contributions to the total pen- 
etration depth, A, when the sample is in the mixed 
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FIG. 5. (Color online) Field versus temperature diagram the 
upper critical field for two different crystal orientations as 
described in the text and also for the location of the maximum 
in x^r (see Fig. 
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FIG. 6. (Color online) Temperature scans of the AC mag- 
netic susceptibility of single crystalline MgCNi3 for various 
values of applied dc magnetic fields up to 8 T. Each run was 
zero-field-cooled and then field-warmed and field-cooled twice, 
which leads to irreversible and reversible portions of the curve. 
Inset: The 6 T temperature scan is used here to define the 
initial irreversible portion of the curve Xirr, the reversible por- 
tion obtained after field- warming and cooling Xrev, and the 
difference Ax at a temperature of 1.4 K. 



state. One of these is the usual London penetration 
depth due to the diamagnetic screening of the applied 
magnetic field by the condensate, XLondon- The other 
component arises from the motion of the vortices and 
a comprehensive expression has been derived in various 
works for \ vortex W^^\t has been shown that in the limit 
of low temperatures and fields that X vor tex reduces to 
the Campbell penetration depthj^l where A^ ompbeH = 
4>oH/a. Taking both contributions into account gives 
the total magnetic penetration depth to be 

A 2 = ^Campbell + ^London- (6) 

In the vortex state, the Campbell penetration depth is 
the dominant term. Here, B is the applied magnetic field 
and </>o is the flux quantum. The Campbell penetration 
depth is important because it contains the necessary in- 
formation to obtain the Labusch parameter, a, which is 
a measure of the curvature of the potential energy asso- 
ciated with the pinning of vortices. 

The dependence of the susceptibility features Xirr, Xrev 
and Ax = Xirr ~ Xrev (defined in Fig. 6) on applied mag- 
netic field can be seen in Fig. 7. It is interesting to note 
the non-monotonic behavior, consistent with a peak in 
Ax near 5 T. This feature, i.e. a maximum in the amount 
of diamagnetic screening at a particular location in the 
H—T phase diagram, is commonly referred to as the peak 
effeclpS and has been observed in other superconductors 
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FIG. 7. (Color online) Upper panel: Susceptibility difference 
Ax = Xirr — Xrev, which can be seen in Fig. 6, at 1.4 K as 
a function of applied dc magnetic field. Lower panel: The 
individual curves of Xirr and Xrev taken at 1.4 K as functions 
of applied dc magnetic field. 



like Lu 2 Fe 3 Si 5 P MgBaP the high-T c BSCCO 2212,^ 
and more recently in iron-pnictides.^^1! Experiments in 
which the resistivity of clean samples of MgCNi3 with 
weak pinning have also shown evidence for the existence 
of the peak effect 

Many explanations have been put forth with the intent 
of explaining the presence of this maximum feature in 
Ax(T, H). One early theoretical work done on the effect 
of disorder induced pinning on a vortex lattice considered 
that instead of the usual Abrikosov lattice, there exists 
a quasi-ordered Bragg glass phase and that the peak ef- 
fect is a sign of the transition from this phase into a 
disordered vortex phaseP^ TDR experiments on BSCCO 
221223 ma y suggest that the observed hysteresis is a re- 
sult of ramping the magnetic field after zero-field cool- 
ing giving rise to macroscopic screening supercurrents, 
j, which shift the vortices into a state of inhomogeneous 
distribution, which is in agreement with the critical state 
(Bean) model. In this scenario, this procedure gives rise 
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to a state consisting of a displaced vortex lattice, which 
disappears when the sample is field-cooled due to a re- 
laxation of screening currents. 

In conclusion, good quality single crystals of MgCNia 
were grown at high pressure and studied using DC and 
AC magnetization. The zero-field London penetration 
depth has been measured and converted into the super- 
fluid density, p s = (A(0)/A(T)) 2 . The conventional weak 
- coupling s— wave BCS temperature dependence of the 
London penetration depth at low temperatures and of 
p s (T) in the whole temperature range can be reproduced 
very well from the TDR measurements with a corre- 
sponding value of A(0) = 245 nm. The H — T phase 
diagram has been mapped by measuring M(T) in differ- 
ent applied dc magnetic fields. The -ff c2 is found to be 
isotropic for two different directions of applied magnetic 
field with -ff C 2(0) « 12.3 T by using the standard Helfand 
and Werthamer analysis. This value corresponds to the 
coherence length of 5.2 nm and together with A(0) = 245 
nm gives a Ginsburg-Landau parameter of k ~ 47. By 



studying the effect of field-cooling versus field-warming 
on the susceptibility, a hysteretic response has been ob- 
served and it has been speculated that this arises due a 
vortex lattice-related phenomenon known commonly as 
the peak effect and signal non - parabolic nature of the 
pinning potential. 



ACKNOWLEDGEMENTS 

We thank V. G. Kogan for useful discussions and H. 
Kim for help with the data analysis. The work at the 
Ames Laboratory was supported by the U.S. Department 
of Energy, Office of Basic Energy Sciences, Division of 
Materials Sciences and Enginnering under contract No. 
DE-AC02-07CH11358. The work at ETH Zurich was 
supported by Swiss National Science Foundation, the Na- 
tional Center of Competence in Research MaNEP (Ma- 
terials with Novel Electronic Properties). 



hxh| 



* |rt-gordon@wi uledu| 
t zhisadlo@phys.ethz.ch 

* wstephen@physik.uz 
§ katrych@phys.ethz.ch 
^ corresponding author: |prozorov@ameslab.gov| 

** work done while at The Ames Laboratory and Department 
of Physics & Astronomy, Iowa State University, Ames, IA 
50011, USA 

" work done while at Laboratory for Solid State Physics, 
ETH Zurich, CH-8093 Zurich, Switzerland 

1 T. He, Q. Huang, A. P. Ramiraz, Y. Wang, K. A. Regan, 
N. Rogado, M. A. Hayward, M. K. Haas, J. J. Slusky, K. 
Inumara, H. W. Zandbergen, N. P. Ong, and R. J. Cava, 
Nature 411, 54 (2001). 

2 J. H. Shim, S. K. Kwon, and B. I. Min, Phys. Rev. B 64, 
180510(R) (2001). 

3 D. J. Singh, and I. I. Mazin, Phys. Rev. B 64, 140507(R) 
(2001). 

4 J. H. Kim, J. S. Ahn, J. Kim, M. S. Park, S. I. Lee, E. J. 
Choi, and S. J. Oh, Phys. Rev. B 66, 172507 (2002). 

5 P. M. Singer, T. Imai, T. He, M. A. Hayward, and R. J. 
Cava, Phys. Rev. Lett 87, 257601 (2001). 

6 L. Shan, H. J. Tao, H. Gao, Z. Z. Li, Z. A Ren, G. C. Che, 
and H. H. Wen, Phys. Rev. B 68, 144510 (2003). 

7 Z. Q. Mao, M. M. Rosario, K. D. Nelson, K. Wu, I. G. 
Deac, P. Schiffer, Y. Liu, T. He, K. A. Regan, and R. J. 
Cava, et al, Phys. Rev. B 67, 094502 (2003). 

8 T. Klimczuk and R. J. Cava, Phys. Rev. B 70, 212514 
(2004). 

9 G. J. MacDougall, R. J. Cava, S. J. Kim, P. L. Russo, A. 
T. Savici, C. R. Wiebe, A. Winkels, Y. J. Uemura, and G. 
M. Luke, Physica B 374-375, 263 (2006). 

10 J. Y. Lin, P. L. Ho, H. L. Huang, P. H. Lin, Y. L. Zhang, 
R. C. Yu, C. Q. Jin, and H. D. Yang, Phys. Rev. B 67, 
052501 (2003). 

11 L. Shan, K. Xia, Z. Y. Liu, H. H. Wen, Z. A. Ren, G. C. 
Che, and Z. X. Zhao, Phys. Rev. B 68, 024523 (2003). 



12 H. S. Lee, D. J. Jang, H. G. Lee, W. Kang, M. H. Cho, and 
S. I. Lee, J. Phys.: Condens. Matter 20, 255222 (2008). 

13 R. Prozorov, A. Snezhko, T. He, and R. Cava, Phys. Rev. 
B 68, 180502(R) (2003). 

14 S. Mollah, J. Phys.: Condens. Matter 16, R1237 (2004). 

15 H.-S. Lee, D.-J. Jang, H.-G. Lee, S.-I. Lee, S.-M. Choi, and 
C.-J. Kim, Adv. Mater. 19, 1807 (2007). 

16 Z. Pribulova, J. Kacmarcfk, C. Marcenat, P. Szabo, T. 
Klein, A. Demuer, P. Rodiere, D. J. Jang, H.-S. Lee, S.-I. 
Lee, and P. Samuely, Phys. Rev. B 83, 104511 (2011). 

17 P. Diener, P. Rodiere, T. Klein, C. Marcenat, J. 
Kacmarcfk, Z. Pribulova, D. J. Jang, H.-S. Lee, H.-G. Lee, 
and S.-I. Lee, Phys. Rev. B 79, 220508(R) (2009). 

18 D.-J. Jang, H.-S. Lee, H.-G. Lee, M.-H. Cho, and S.-I. Lee, 
Phys. Rev. Lett. 103, 047003 (2009). 

19 T. G. Amos, Q. Huang, J. W. Lynn, T. He, and R. J. Cava, 
Solid State Commun. 121 73 (2002). 

20 G. Sheldrick, SHELXS-97, Program for the Solution of 
Crystal Structures (University of Gottingen, Germany, 
1997); SHELXL-97, Program for the Refinement of Crystal 
Structures (University of Gottingen, Germany, 1997). 

21 N. D. Zhigadlo, S. Katrych, M. Bendele, P. J. W. Moll, 
M. Tortello, S. Weyeneth, V. Yu. Pomjakushin, J. Kanter, 
R. Puzniak, Z. Bukowski, H. Keller, R. S. Gonnelli, R. 
Khasanov, J. Karpinski, and B. Batlogg, Phys. Rev. B 84, 
134526 (2011). 

22 A. Kumar, R. Jha, S. K. Singh, J. Kumar, P. K. Ahluwaha, 
R. P. Tandon, and V. P. S. Awana, J. Appl. Phys. 11, 
033907 (2012). 

23 R. Prozorov and R. W. Giannetta, Superc. Sci. Technol. 
19, R41 (2006). 

24 R. Prozorov and V. G. Kogan, Rep. Prog. Phys. 74, 124505 
(2011). 

25 R. Prozorov, R. W. Giannetta, A. Carrington, and F. M. 
Araujo-Moreira, Phys. Rev. B 62, 115 (2000). 

26 H. Ki m, N. H. Su ng, B. K. Cho, M. A. Tanatar, R. Pro- 
zorov, |arXiv:1301.6"054"l (2013). 



8 



27 E. Helfand and N. R. Werthamer, Phys. Rev. 147, 288 
(1966). 

28 R. Prozorov, R. W. Giannetta, N. Kameda, T. Tamegai, 
J. A. Schlueter, and P. Fournier, Phys. Rev. B 67, 184501 
(2003). 

29 M. W. Coffey and J. R. Clem, Phys. Rev. Lett. 67 386 
(1991). 

30 E. H. Brandt, Phys. Rev. Lett. 67, 2219 (1991). 

31 C. J. van der Beek, V. B. Geshkenbein, and V. M. Vinokur, 
Phys. Rev. B 48, 3393 (1993). 

32 A. M. Campbell, J. Phys. C, 2 1492 (1969). 

33 W. DeSorbo, Rev. Mod. Phys. 36, 90 (1964). 

34 R. T. Gordon, M. D. Vannette, C. Martin, Y. Nakajima, 
T. Tamegai, and R. Prozorov, Phys. Rev. B 78, 024514 



(2008). 

C. Martin, M. D. Vannette, R. T. Gordon, R. Prozorov, N. 

D. Zhigadlo, and J. Karpinski, Phys. Rev. B 78, 144512 
(2008). 

P. J. W. Moll, R. Puzniak, F. Balakirev, K. Rogacki, J. 
Karpinski, N. D. Zhigadlo, and B. Batlogg, Nature Mate- 
rials 9, 628 (2010). 

N. D. Zhigadlo, S. Katrych, S. Weyeneth, R. Puzniak, P. 

J. W. Moll, Z. Bukowski, J. Karpinski, H. Keller, and B. 

Batlogg, Phys. Rev. B 82, 064517 (2010). 

P. Prommapan, M. A. Tanatar, B. Lee, S. Khim, K. H. 

Kim, and R. Prozorov, Phys. Rev. B84, 060509 (2011). 

T. Giamarchi and P. Le Doussal, Phys. Rev. B 52, 1242 

(1995). 



